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Abstract 

The isotropic S = 1/2 antiferromagnetic spin ladder with additional four-spin ring exchange is 
studied perturbatively in the strong coupling regime with the help of cluster expansion technique, 
and by means of bosonization in the weak coupling limit. It is found that a sufficiently large 
strength of ring exchange leads to a second-order phase transition, and the shape of the boundary 
in the vicinity of the known exact transition point is obtained. The critical exponent for the gap 
is found to be n — 1, in agreement both with exact results available for the dimer line and with 
the bosonization analysis. The phase we determined after crossing the critical line is gapped and 
spontaneously dimerized. The results for the transition line from strong coupling and from weak 
coupling match with each other naturally. 

PACS numbers: 75.10.Jm, 75.40. Cx, 75.40.Gb, 75.30.Kz 



1 



I. INTRODUCTION 



At half filling and in the limit of small ratio x = t/U of hopping and on-site Coulomb 
repulsion the Hubbard model can be mapped to an effective spin exchange Hamiltonian. In 
the leading order in x the standard (bilinear) antiferromagnetic nearest-neighbour Heisen- 
berg exchange interaction with the exchange constant J = 2t 2 /U is obtained. Terms of 
higher order in x yield, except bilinear exchange interactions beyond the nearest neighbors, 
also exchange terms containing a product of four or more spin operatorsBSi'S. Those higher- 
order terms were routinely neglected up to recent time, when it was realized that they can 
be important for a correct description of many physical systems. 

For the first time biquadratic exchange was used for the description of the magnetic 
properties of solid 3 Hel. Recently it was suggested that some strongly correlated electron 
systems like cupratesi'l and spin ladders!'! are expected to exhibit ring exchange. The 
analysis of the low-lying excitation spectrum of the p-d-model shows that the Hamiltonian 
describing Cu02 planes should contain a finite value of ring exchangeSJli'00. 

There is a number of experimental work like inelastic neutron scatteringi!! and nuclear 
magnetic resonanc ,00 on Sr 14 Cu 24 04i and Ca 8 La 6 Cu240 4 i as well as optical conductivity 
measurement s0'0 on (Ca,La)i4Cii2404i. All these substances contain spin ladders built of 
Cu atoms. The attempts to fit the experimental data without taking ring exchange into 
account yield an unnaturally large ratiolii of Jw/Xune ~ 2 which is expected neither from 
the g eo m etnca, structure of the .adder nor froufrou, structure ca.cu.atioofl It can 
be shown that inclusion of other types of interactions in particular an additional diagonal 
interaction does not remove this discrepancy!. 

In the present paper we study the isotropic 5 = 1/2 antiferromagnetic spin ladder with 
additional four-spin ring exchange. Starting from the dimer limit of uncoupled rungs, we 
use the cluster expansion technique to calculate the dispersion of the elementary excitations 
to high order in the perturbation parameters J\ eg and J r i ng , and study the influence of the 
four-spin exchange on the spectrum. In the (Ji eg , ./ring) space we have found a transition 
boundary where the gap vanishes. The shape of this boundary is obtained in the vicinity 
of the exactly known transition point lying on the "dimer line" where the exact ground 
state is a product of dimersll. With the help of Pade-approximants we calculate the critical 
exponent rj for the gap and obtain r\ ~ 1, which is in agreement with the exact results 
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available for the dimer line. We also analyze the opposite limit of weakly coupled chains by 
means of the bosonization technique and come to the same conclusion on the linear behavior 
of the gap. By this method we investigate the other side of the critical line. We show that 
the phase which emerges above the critical value of J r i ng is a spontaneously dimerized phase 
with a finite gap to the elementary excitations. This result is supported by recent numerical 
calculations 



II. MODEL 



We consider the isotropic S = 1/2 antiferromagnetic spin ladder with additional cyclic 
four spin exchange. Fig. [I] illustrates the Hamiltonian, which is of the form 



Tt — 7i TUllg + Ti.\ cg + Ji 
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In (JTJ) (ijkl) labels a four spin plaquette. Piju leads to a cyclic permutation of spin moments, 
i.e. 
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We rewrite the operator Pijki as a product of two spin permutation operators and obtain the 
following result which contains both bilinear and biquadratic terms of the spin-l/2-operators: 
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In further discussion the constant term is omitted, and periodic boundary conditions are 
used. 

On the so called "dimer line" where J ring = Ji eg the ground state of the model (|1|) is a 
product of singlet dimers placed on the ladder rungsJUll Moreover, on this dimer line the 
propagating triplet becomes an exact excitation,!^ and its energy is given by the following 
simple expression (the lattice constant is set to unity): 



£t(q) = J 



— 11- 

rung ^ ° ring 



2 Jring COs(g) 



(4) 



It is easy to see that the gap of this excitation vanishes at q = 7i for J ring = | J rung . It was 
suggested,i'0 that this point belongs to the line which has to be identified as corresponding to 
the transition into the dimerized phase with spontaneously broken translational symmetry 
along the ladder. However, recent numerical calculations^!! have created some doubt, 
indicating the possible existence of a gapless phase on the other side of the transition line. 

Si,;-i j lpo . Si |( j Si jt _|_i Sx,»+2 




>2,i+2 



FIG. 1: Schematic structure of a two-leg ladder with additional ring exchange 



III. ELEMENTARY EXCITATIONS AND PHASE DIAGRAM 
A. Rung dimer limit: Cluster expansion 

In this section we study the low-lying excitations of the above model by perturbation 
theory. Therefore we start in the dimer limit where J ruri g is the only non- vanishing exchange 
constant and measure all other interactions in units of J Tnng . In fact we have only two 
coupling parameters: 

Q?L Jleg/ Jrung and Gift Jring/ Jrung- (5) 
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The dispersion is a function of both a R and «l possessing following form: 



oo n n 



^(g)/ Jmng = Yl Yl C< Cj a L a R 3 COS (^)- ( 6 ) 

n=0 i=0 j=0 

The dispersion was obtained by using the cluster expansion formulationHEil of the per- 
turbation theory in a R , oil- Clusters with maximum of 14 edges were considered, which is 
an improvement of 11 orders to our previous work!. Up to the third order one has 



cu (3) (g) = fi + /ii cos(g) + fi 2 cos(2g) + /i 3 cos(3g) 

3 

Ho = 1 - 2ct R + - (a L - a R ) 2 (2 + a R + 5 a L ) 



A*i = «r + «l - ^ («l - or) 2 («r + «l) (7) 

1 2 

yu 2 = -- («l - a R ) (1 + a R + «l) 
A^3 = ^ ("l - or) 2 («r + a L ) • 

o 

One may notice that for a-^ = a R only the first-order term is left in the expansion, and one 
recovers the exact result (|J). 

In Fig. and Fig. [3] we plot some typical dispersions for fixed «l = 0.25 («l = 0.50) and 
varying a R . One can see that the gap at q = it decreases with increasing a R and finally 
vanishes. It is also seen that the convergence of the method becomes worse when «l and 
a R increase. However, the results from direct series expansion can be improved by using the 
Pade- approximation technique. 

We have studied the vicinity of the exactly known transition point «l = «r = 0.25 
and calculated the phase boundary where the gap closes. The resulting view of the phase 
diagram is shown in Fig. |5](a). We were not able to proceed beyond the intervals which are 
marked by the arrows by increasing the order. It seems that we hit the convergence radius 
of the present method. 

We have also calculated the critical exponent r\ of the gap A ring as a function of both 
and a R where we use the following definition: 

A ring oc \a L - a c L \ v{ai{) and (8) 
A ring ex \a R -a R \^\ (9) 

Within the convergence interval «l = [0.17,0.36] resp. a R = [0.24,0.38] this exponent is 
equal to 1 with the accuracy of ±0.01. This is in agreement with the exact result (HI) and 



shows that the picture of the phase transition does not change when one moves away from 
the exact point «l = «r = 0.25. This behavior of the gap indicates that the phase on 
the other side of the phase boundary is gapped. Indeed, linear growth of the gap indicates 
presence of a relevant operator, whose amplitude changes sign at the boundary. Existence 
of the high-«R gapless phase would, in contrast, imply presence of a marginal operator 
which becomes irrelevant on the boundary. That would necessarily mean exponentially slow 
growth of the gap, which contradicts to our observations. Those numerical results are also 
in a good agreement with the conclusions of bosonization analysis (see the next section), 
which predict a linear behavior of the gap at both sides of the transition. 



B. Decoupled chains limit: Bosonization 

In this section we perform the weak-coupling analysis (J\ cg ^> J r i ng , J rU ng) of our model. 
In a weak coupling-bosonization analysis we find a critical line separating rung singlet phase 
from spontaneously dimerized phase. We conclude that the critical line as well as the phases 
it separates match with each other in weak and strong coupling regimes. We show that the 
effect of the ring exchange in the infrared limit is reduced (up to marginal corrections) to 
that of the leg-leg four-spin interaction as considered earlier by Nersesyan and TsvelilH m 
the context of spin-phonon coupling. 

We write the spin operators on each chain in terms of their smooth and staggered parts: 

Si, 2 (x) = J 1)2 {x) + {-l) x n 1)2 {x). (10) 

The effect of rung interaction is well understood by bosonization using Majorana fermion 
formalisml!. Now we will bosonize ring exchange, treating it as a perturbation to decou- 
pled chains and at the end we will add bosonized-refermionized terms coming from rung 
interaction as in@. We decompose the Hamiltonian in the following way: 

H(x) = H^x) + H 2 {x) + H quad {x) + H hiqnad {x) , (11) 

where Hi 2 are critical Gaussian models describing first and second decoupled chains. H qnad 
stands for the quadratic spin interactions and -ffbiquad for the four-spin interactions originat- 
ing from the ring exchange term. 
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Neglecting renormalization of the intrachain interaction we first analyze the quadratic 
spin interactions which can be cast in the following form: 

-Hquad ~ (^ring + ^ring 

)Ji(x)J 2 (x) 

+ (^ng-^rin g )ni(a:)n 2 (x) (12) 

where J r f ng = J r * = J T i ng . Since the scaling dimension of the smooth part of the spin oper- 
ator is 1, while the dimension of the staggered part is 1/2, no relevant terms are generated 
from the quadratic part of the ring exchange, and only marginal terms are left. After bosoniz- 
ing the biquadratic part only leg-leg interaction will survive, because diagonal-diagonal and 
rung-rung terms give non-distinguishable relevant contributions in the infrared limit which 
cancel each other due to the overall opposite signs in front of them (which is fixed by the 
structure of the ring exchange): 

T^biquad ~ (^ring ^ring ^ring 

)e 1 (x)e 2 {x), (13) 

where 61,2(2) = (— l^Si^a^Si^x + ao) represents the dimerization operator of the first and 
second chain, respectively, and = J^ g = J^ g = 2 J ring . In the Majorana representation 

(retaining only relevant operators) we arrive at the following Hamiltonian: 

3 

/ — in 

(14) 



H = ^T I dx 



a=0 

.2, 
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with m = — c 2 a J r i n g/27r. c is the Lukyanov-Zamolodchikov constant at the 577(2) AFM point 
and a is a non- universal, cutoff-dependent positive constant. Thus in the weak-coupling 
limit we have effectively reduced the ring exchange to the leg-leg interaction. The only 
difference between the bosonized forms of the ring exchange and the pure leg-leg coupling 
stems from the marginal current- current interaction which does not appear in the leg-leg 
biquadratic interaction. In contrast to the leg-leg interaction ring exchange is not invariant 
under independent global SU(2) rotations of spins on each chain, and thus should not 
enjoy full 0(4) symmetry. This symmetry is in fact lowered by marginal operators. The 
refermionized version of the marginal current-current interaction contained in ring exchange 
will take the following form in the Majorana representation:!^ 

H ma , rg = <7ring«0 J ^[(£r6z,) (CflCi) 

+ (elei)(ele!) + (eki)(elel) 

- (e R e L +e R e L +e R e L mi)] as) 



Renormalizing the masses it weakly splits the 0(4) quadruplet into a triplet and a singlet, 
consistent with the symmetries of ring exchange: 

m t Jring&o . Jieg 

— — > 1H — l n i — 7) 

m ttv \m\ 

m s y ^ _ 3^ngO0 ^lcg (\Q) 

m ttv \m\ 

The full refermionized model, including both rung and ring exchange, in the Majorana 
representation reads as: 

r -ivt , 



H = dx 



a=l,2,3 

dx 



where 



c 2 

IT^t Q ("^rung ^"/ring)) 
C 2 



(17) 



m s = - — (3J rung + aJ r i ng ) (18) 
From the above formulae, one readily obtains the line where the triplet mass vanishes. Ac- 



cording to Ref. [25|, on this fine-tuned line a phase transition from conventional Haldane 
phase (rung exchange dominated phase) where spectrum displays coherent single-particle 
(magnon) excitations to non-Haldane spontaneously dimerized phase (ring exchange dom- 
inated phase) without coherent magnon modes takes place. This transition belongs to the 
universality class of critical, exactly integrable, S = 1 spin chain (Takhtajan-Babujian point) 
with the central charge c = 3/2. The dimerization pattern emerging after crossing the crit- 
ical line is the following: the chains become dimerized in a staggered way to each other 
with a nonzero relative dimerization. This is consistent with the fact that for the anti- 
ferromagnetic interchain interaction effective S = 1 spins exhibiting nonzero string order 
are formed across the ladder diagonals rather than along the rungs. On either side of this 
line the system is gapped, described in terms of free massive Majorana fermions with the 
symmetry SU(2) <g> Z 2 . The gap (which is the mass of the Majorana triplets) opens linearly 
as one deviates from the criticality. Owing to the £77(2) symmetry of the model no other 
perturbations than mass terms of Majoranas are allowed. In the weak-coupling limit the 
existence of gapless excitations other than on this line is thus excluded. 
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We argue that this quantum critical line (as well as the phases it separates) determined 
in the weak-coupling limit is smoothly connected to the one discovered in the strong rung 
coupling limit. The existence of the lowest triplet excitation gap which vanishes linearly as 
we approach the critical line from the rung dimerized side of the phase diagram (which is 
seen in both the weak coupling and the strong coupling limit) is taken as evidence that this 
part as well as the critical line are smoothly connected in weak and strong coupling limit. 

A closer inspection of our results shows that the smooth connection of the transition line 
in the two regimes is verified quantitatively: the numerical data shown in Fig. [|(a) seem to 
indicate asymptotic behavior of the phase boundary in the limit of large J\ eg , 

lim 7 C IP m 22 

This number is consistent with the data for the critical line obtained from Lanczos exact 
diagonalizationi for Ji og /Jmng ~ 1 arid for finite system size (N = 24). In an alternative 
presentation of J r i ng vs. J rU ng (see Fig. |(b)) this asymptotic behavior translates into a 
straight line joining smoothly the result following (|TTD, J r c ing ~ J r c un g- Thus the limiting 
value J^ing/ ^rung ~ 0-^2 as obtained in strong coupling appears at the same time a numerical 
determination of a -1 in 

Based on that the critical line is specified in terms of three massless Majoranas: the 
SU(2)2 Wess-Zumino model is the universality class describing this line throughout the 
phase diagram (consistent with numerical results^). The gapless behavior arises due to 
fine tuning of rung and ring exchange parameters which leads to mass cancellation of the 
Majoranas. Upsetting of this fine tuning on either side of the critical line will give rise to 
nonzero Majorana masses. 



IV. CONCLUSION 



In summary, we have studied the isotropic S = 1/2 antiferromagnetic spin ladder with 
four-spin ring exchange, both in the limit of weakly coupled chains and in the strongly 
coupled (dimer) limit. 

In the dimer limit, the use of the linked cluster expansion technique has allowed us to 
calculate the dispersion of the elementary excitations to a high order in the perturbation 
parameters J\ eg and J r i ng , and to establish the shape of the transition line where the gap 
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vanishes linearly. In the opposite limit of weakly coupled chains, using the continuum-limit 
bosonized form of the Hamiltonian, we have come to the conclusion that the gap behaves 
linearly at both sides of the transition. We identify the phase which emerges above the 
critical value of J r i ng as a spontaneously dimerized "non-Haldane" phaseil whose elementary 
excitations are pairs of massive kinks. From our two complementary perturbation approaches 
a consistent description of the unique transition line in the full phase diagram has emerged. 
However our methods cannot access the regions where J mng ~ J\ eg and J r i ng > J rU n g , Jieg- 
Therefore it is not excluded that some other phases emerge within those regions as indicated 
by DMRG calculations^. 
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FIG. 2: Dispersion for = 0.25 up to the 14 th order and for various oil- On the exact line 
(«r = 0.25) the gap vanishes. Further increase of or > 0.25 leads to a loss of convergence. 
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FIG. 3: Dispersion for «l = 0.50 up to the 14 th order and for various «l- 
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FIG. 4: The critical exponent 77 as a function of q?l resp. or using a log 2 [6, 6] Pade approximant. 
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FIG. 5: Phase diagrams: (a) result from the 14 th order of perturbation theory in or = Jrmg/^mng 
and «l = ^leg/^rung- The arrows mark the regions where the perturbation theory is considered to 
be valid, while the dotteded and dashed parts of the phase boundary indicate the regions where 
the cluster method loses convergence; (b) the SclIXLG 111 ( taring i^rung) P&I*ciXfl6t6r 

space showing the 

smooth connection between bosonization results at small and numerical results at large J r i ng and 
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